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Abstract 

We have developed new methods to calculate dispersion curves (analytically 
in the simpler cases) from which we are able to derive the spatial distribution 
of electron and current densities. We investigate the case where the magnetic 
field varies linearly with position and the results provide useful insights into 
the properties of this and other field distributions. We consider spin as well as 
a confining electrostatic potential. We show that the electron and the current 
density exhibit a very rich structure related to the quantisation of the energy. 
Moreover there is a direct contribution to the current density due to the spin 
which could be of interest in relation to spin polarised current. 

PACS numbers: 73.20.Dx, 73.40.-c, 73.40.Hm 
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I. INTRODUCTION 



A 2-dimensional electron gas (2DEG) in a magnetic field has proved to be an extremely 
rich subject for theoretical and experimental investigation For example, considerable 
effort has been devoted to the study of the Integral and Fractional Quantum Hall Effects (I 
and FQHE), transport properties and edge states 0. Except in a few cases however, 
the magnetic field considered was homogeneous. In this article we address the problem of 
a magnetic field varying linearly with position added to an underlying homogeneous field. 
This is of relevance because in real systems: (a) a constant magnetic field is not always 
attainable, and (b) an inhomogeneous field may be desired. Another interesting point is 
that in the Composite Fermion (CF) theory |^,^ which is used to describe the FQHE, the 
electron-electron interaction, necessary for the appearance of the FQHE, is incorporated 
into an effective magnetic field via a singular gauge transformation. The result is a system 
of non-interacting quasi-particles carrying a fictitious magnetic flux in an inhomogeneous 
effective magnetic field. A better understanding of the properties of a simple non-interacting 
electron gas in an inhomogeneous field might therefore bring useful insights into CF theory. 



II. MODEL 

To investigate the electronic properties of a non-interacting electron gas in a linearly 
varying magnetic field we consider the following Hamiltonian 

where A = (J^Biy'^ + B^y, 0, 0), B = —(0, 0, Biy + Bq), the confining potential due to the 
walls Vc{y) = (3 exp (^(y — 1/e)) + exp + Z/e)) , Ve is the position of the edge of the 

system, S is the spin operator and g = 2 for the g-factor. The parameters a and (3 allow the 
shape of the potential to change continuously from very sharp to very smooth, which can 
modify the properties of the system MM. By solving 
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we can obtain the electron density 



p{x,y) = x*{x,y)x{x,y), 



(3) 



states 



(4) 



where we sum over all states with energy E < Ep, and the current density for a state n 

y) = -^Re [Xn{x, y) (p - eA) Xn{.x, y)] + —V x [x*n{x, y)Sxn{x, y)] . 

As a result of our choice of A and B, the symmetry of the Hamiltonian allows us to write 
the wave function as x{x,y) = e'^'^'^^ip{y), where kx is a good quantum number, and we then 
obtain the 1-dimensional Schrodinger equation 



. (hkx - eAxf - —B, + Vc 
2m* 2m* 



i){y) = Eip{y), 



(5) 



with s = ±0.5. Eq.(5) then enables us to derive the dispersion curves En{kx), and to rewrite 
the electron density as 



Piy) = Y.Erniy)My) 



(6) 



n 



and the current density for one state n and a fixed kx as 



(7) 



me 

both now solely functions of y. Due to the symmetry of the system there is a current density 
only along the x axis. Integrating over y will give the total current /" carried, for a fixed kx 
by the state n and summing over n and kx the total current Ix- 



III. METHOD 



Starting from Eq.(5) the Hamiltonian can be written as 



Py m*uji A m*uJiUJo q / m*ujn hkxOJi \ ^ 
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+ Ve{y) 



ip{y) = Eip{y), 



(8) 
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with luq 



eBn 



^ and coi = For what follows it is useful to introduce the dimensionless 



variable y — (h ^m*Oi)ij ^ y and p — (rn*ujih^^ ^ p which yields for (8) 

\^{y)^E^{y), 



^ + aj/^ + by^ + cf + dy + e + Kl 



(9) 



with m' 



(m*a)^fi)V3 ^1 ^' ^ given in unit of energy. Although in some simplified 

cases it is possible to obtain analytical results, as we will see below, there is in general no 
way to find the analytical solution of Eq. (9), and therefore we have to resort to numerical 
calculations. Eq.(9) can be solved by expanding in terms of oscillator functions, 0n(y) — 
H„(y)e 2 where H„ is a Hermite polynomial, and then by numerically diagonalising the 
corresponding secular equation 



Detliffen - E5kn\ = 0. 



(10) 



Using the properties of the Hermite polynomials all the matrix elements can 
be calculated analytically (Appendix) which greatly improves the diagonalisation method. 
However before starting with the numerical calculations we can try an analytical approach 
to Eq.(8) in the simplified case where = and Vc — 0. We then have 
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2m* ■ 2m* V 4 h"^ ^ % ' km. 
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We choose two regimes for which Bi ^ 0: (a) < 0, single well potential (SWP), near 



We 



S = 0, and (b) k^: > 0, double well potential (DWP), near B = ±Bi ^j2hk^/eBi 
expand parabolically around the minima of the effective potential and obtain harmonic 
oscillator equations which we solve analytically. The expressions obtained for the energy are 



SWP : Er, 



2m* 



+ (2n+ l)i 



'kxsBi eBiT] 



n 



DWP : Er, 



2m* 



(2n+l)J^-^±2,./2^^^^^ 



(12) 



(13) 



h 2hkx V h 

where rj = From here it is straightforward to derive the group velocity (^^f^) for the 
state n as 
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SWP : Vr 



kx + 



eBiff 
2hkl 



(2n + l)i 



' eBi 
IQhkr 



(14) 



DWP : Vr 



8hkx 



2hk 



(15) 



IV. RESULTS 

In the following calculations the system we consider corresponds to an ideal slab of 
GaAs/AlxGai_xAs heterostructure filled with an ideal 2DEG. The effective electron mass is 
m* = 0.067 me, the electron density 4 x 10~^j4 ^, and the sample has width, when 7^ 0, 
2?/e = 2 X 10^ j4 and length L ^ 1 with periodic boundary conditions along x. 

In Fig.|^ are reported the numerical and the analytical results (Eqs.(12),(13)) for the case 
when Vc = 0, Bi = IG/A, and Bq = 0. We can see that the agreement is very good except 
around k^ = where the method breaks down. Although it is possible to obtain useful 
information from analytical calculations they do not allow us to derive complete dispersion 
curves and hence the electron or current densities. Moreover we are interested in taking into 
account the effects of a confining potential Vc but this cannot be included in our analytical 
approach. We have then to use a numerical approach. 

So, starting from the same model system configuration as above but this time with 
Vc ^ and, Bq = 0.5 * 10^ G. We chose a = 100 and (3 = 50 which correspond to quite a 
sharp confining potential. The dispersion curves are plotted in Fig. 0. The degeneracy of 
the energy levels is completely removed and the new structure appearing in the dispersion 
curves is due to the breaking of the y symmetry in the Hamiltonian. It is interesting to 
note that some similar features were observed in the case of a curved 2DEG in a constant 
magnetic field 0. Having En and ipn we can now calculate j^\y) and p{y), but to do this 
we need the Fermi energy Ep. This can be obtained by minimising the total energy with 
the constraint that the number of electrons is constant with A^ given by 

N = — yf;'dk^^\ (16) 
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where {k^ } are the parameters we vary to minimize the energy. As we might expect, Ep 
does not depend on but, in contrast to the assumption in 0, it is not independent of the 
magnetic field. This is shown in Fig. ^ where Ep has been calculated for various values of 
Bi and Bq. This dependence is due, in the absence of external leads, to the walls, which can 
be seen by the fact that when Bi, which in contrast to Bq removes the degeneracy of the 
states and gives rise to an effective confining potential, increases, Ep becomes independent 
of the magnetic field. 

Using Ep we can now calculate the electron density p{y)- In Fig. § p{y) is plotted for 
different values of Bi and Bq. We see that with an external confining potential (V^ 7^ 0) 
p{y) is not constant; it can have quite a rich structure with local charging effects. Moreover, 
when Bq 7^ 0, p{y) becomes asymmetric in y. It has to be noted that in the case K = p{y) 
is constant the oscillations on left and right hand side of the graphic are only a numerical 
effect due to the fact that in this case one should use a larger set of (pn for the expansion of 
ip- But for Vc 7^ the structure of p{y) can be explained in the following way. When the 
number of states as a function of k^ is plotted a discontinuous curve due to the quantisation 
of the energy is obtained. Although there is no simple relation between the k^ space and the 
real space, as for a homogeneous magnetic field, the number of oscillations of the electron 
density is the same as the number of steps in the discontinuous curve and then is indeed a 
consequence of the energy quantisation. 

This can be seen in Figs. |^ and The next step now is, by means of Eq.(7), to calculate 
the total current density given by 

jM = E itHy), (17) 

states 

where the sum runs over all the states with E < Ep. In Fig. |^ are reported the current 
densities for different magnetic fields disregarding first the part due to the spin. One sees 
that the shape of the current density is much more subtle than might be expected from 
semi-classical approximations. Actually considering j^\y) for \kx\ ^ or, in other words 
for the largest kx close to the Fermi energy, the movement of the electron can be described by 
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its classical orbits (drifting orbits along the edge and snake orbits in the opposite direction 
along the line where 5^0). For smaller k^, however, the states tunnelling between the 
two wells but mainly the ones of energy above the central maximum of the double well are 
very important and their contribution cannot be overlooked. In fact j^\y)^ after summing 
over all n and A;^, turns out to be very small for the case where -Bq = for U- In Fig- |§ 
are reported different i^\y)- For = — O.OlSyl ^ , k^ = 0.08^4 ^ and n = the movement 
of the electron is well defined by its classical orbits (edge and snake orbits). But around 
kx = and for exemple here n = 31 the situation is more complicated. An interesting point 
is that now the current density for one state can be positive and negative as a function of y. 
Moreover we see that the positive part is located in y where the density of current flows in the 
other direction due to the presence of the snake orbits. The same kind of phenomena appears 
with edge orbits. This can be understand by considering the first term hk^ — \eBiy^ — eB^y 
in Eq.(7). It is easy to imagine that when summing the current density over all the different 
states the result is rather different from what we might expect from the consideration of the 
simple classical picture of the orbits. It has to be noted that although the current density 
can be positive and negative as a function of the group velocity Vg = J dy j^\y) has a 
well defined sign and has been verified from the dispersion curves by means of the relation 
Vg = ^^f^- In order to get a better understanding of the shape of the current density we 
have plotted together the electron and the current densities in Fig. |^ andp!0|. The density 
of current oscillates between positive and negative "channels" as a function of the electron 
density and then is a reflection of the quantisation of the energy. It is interesting to note 
that in the study of the FQHE there is also the appearance of channels which can be seen 
there as alternating strips of compressible and incompressible fluid ||7|JT0|]. When Bq ^ 



the current density increases with Bq and flows in opposite directions on both sides of the 
sample when Bq is large enough compared to Bi or in other words when Bq is large enough 
to overcome the effective confining potential due to Bi so that the electrons are confined by 
the external confining potential Vc. Moreover the current density becomes asymmetric. 
Until now only the first term in Eq.(7) has been considered but there is a second term 



containing the derivative of the electron density and directly related to the presence of 
the spin for the electron. Because, as we have seen above, the electron density displays 
a very rich structure, one can expect some contribution to the current density due to the 



spin of the electron. This is shown in Figs. |Tl| and |1^. Although the part due to the 
spin is smaller than the first term in Eq.(7) it is nevertheless noticable. This could imply 
some interesting phenomena in relation to spin polarised currents. The problem is that the 
"channels" are more or less at the same position for spin up or down which makes it quite 
difficult to distinguish between both spin directions. On the other hand, many parameters 
can be varied, such as the magnetic field, the width of the system or the external confining 
potential which may allow us to find a suitable system configuration for the production of 
polarised currents. 

Finally it has to be stressed that all the discussion above concerned the current density. 
Although this quantity is non-zero and has a rich structure, it does not imply that the net 
current = J2kx 1 3x{y) dy is non-zero. In fact as our calculation showed Ep is independent 
of kx which means there is no difference of potential across our system and thus no net 
current. 

V. CONCLUSIONS 



In this work we have studied the effect of a linear magnetic field in a 2DEG. In certain 
simplified cases we were able to carry out some analytical calculations and to derive the 
dispersion curves En{kx) for quite a large range oik^ ^ These results were found to be in 
very good agreement with our numerical results. In the general case with an external confin- 
ing potential we carried out numerical calculations. We derived the whole of the dispersion 
curve and using it, we calculated the electron and current densities. It is worthwhile noting 
that for this calculation we need to consider the states for /cj, > as well as states for A;^. < 
and kx — 0. This point is important and could help in understanding some recent results |TT 



obtained in the framework of CF theory. For the derivation of the electron and current den- 
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sity we first calculated the Fermi energy Ep taking into account that the number of electrons 
is constant. It turned out that although Ep is independent of it is, however, a function of 
the magnetic field. This is due to the external confining potential. The electron and current 
densities show a very rich structure which can be seen as a consequence of the quantisation 
of the energy, although there is no simple relation between the kr^ and y space as is the case 
for a constant magnetic field. Moreover the current density exhibits alternating "channels" 
of positive and negative current. It would be interesting too to include interaction between 
electrons and to study the effect of the self-consistency on the way the energy levels cross 
the Fermi energy. This can give rise to interesting phenomena, particularly in connection 
with the shape of the electron density [[7|,p!2|. Finally, because p{y) is not constant we have 
a contribution to the current density directly due to the spin of the electron which could 
imply some interesting phenomena in relation to spin polarised currents. 
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APPENDIX A: 



Using the properties of the Hermite polynomial the matrix elements Hkn can be derived 
in a recursive way. {(pkl + ay^ + by^ + cy'^ + dy + \4>n) gives 



n — k 



2m' V2 y V \ 2/4/ V2 



n = k + 1 



n^k + 2 



'{k + i){k + 2y 



n — k + 3 



' {k + l){k + 2){k + 3) y ^ 



k + A 

f{k + l){k + 2){k + 3){k + A)V^ 
16 " 



and (0n|K|0m) due to the confining potential 
n — k 



n — k + i 

^ \n{n-l)...{n-i + l)) \2aJ " [ 2a'' I I ^ ^ ^ 



1 



with the associate Laguerre polynomial. 
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FIGURES 

FIG. 1. Dispersion curves obtained analytically compared with numerical results. Vc = 0, 

Bi = 1 G/ A and Bq = 0. The curves plotted correspond to the levels n = 0, 1, 10, 11. 

FIG. 2. Dispersion curves with ^ 0, 5i = 1 G/A and Bq = 0.5 x 10^ G. 

FIG. 3. Fermi energy as a function of the magnetic field B. 

FIG. 4. Electron densities for different magnetic fields. The solid line corresponds to the case 
without an external confining potential. 

FIG. 5. Electron density p(y) ioi Bi = 1 G / A and = 0. The insert shows the number of 
states as a function of kx- The number of oscillations in p{y) corresponds to the number of states. 

FIG. 6. Same as Fig. 5 but with Bi = 0.5 G/A and Bq = 2x 10^ G. 

FIG. 7. Current density for two different magnetic fields. — : Bi = 1 G / A, Bq = Q and • • • : 

Bi = 0.5 G/A, Bo = 2x 10^ G. 

FIG. 8. Current density calculated for a fixed kx and a fixed n. For kx = —0.018 A ^ (■ ■ •), 
kx = 0.08 A ^ (--), and n = the movement of the electron is well defined by its classical orbits 
(edge and snake orbits). Around kx = 0, here kx = 0.005 A ^ and n = 31 ( — ), the situation is 
more complicated. The density of current can be positive and negative as a function of y. 

FIG. 9. Total density of current (sum over all states with E < Ep) for Bi = 0.5 G/A, 
Bq = 2 X 10^ G. The density of current ( — ) oscillates between positive and negative "channels" 
as a function of the electron density (•••)• 

FIG. 10. Same as Fig. 9 but with Bi = 3 G/A, Bq = 2 x 10^ G. 
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FIG. 11. Contribution due to the spin to the current density. The dotted line is the other 
contribution in Eq.(7) to the current density. Bi = 0.5 G/A, Bq = 2 x 10^ G. 

FIG. 12. Same as Fig. 11 but with Bi = 3 G/A, Bq = 2 x 10^ G. 
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